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Abstract 

In the two-Higgs-doublet model the alignment of the Yukawa matrices in flavour space 
guarantees the absence of flavour-changing neutral currents at tree level, while introduc- 
ing new sources for CP violation parametrized in a very economical way [1] . This implies 
potentially large influence in a number of processes, b — > sj being a prominent example 
where rather high experimental and theoretical precision meet. We analyze the CP rate 
asymmetry in this inclusive decay and determine the resulting constraints on the model 
parameters. We demonstrate the compatibility with previously obtained limits [2]. More- 
over we extend the phenomenological analysis of the branching ratio, and examine the 
influence of resulting correlations on the like-sign dimuon charge asymmetry in B decays. 



1 Introduction 



Adding a second doublet to the scalar sector of the Standard Model (SM) has been a popular 
idea since its proposal in the 70's [3], due to its simplicity and being the low-energy limit 
of some more complete theories. In the most general version of the model, the fermionic 
couplings of the neutral scalars are non-diagonal in flavour and, therefore, generate unwanted 
flavour-changing neutral-current (FCNC) phenomena. Different ways to suppress these have 
been developed, giving rise to a variety of specific implementations [31-119]. A new approach 
in that respect is given by the Aligned Two-Higgs-Doublet Model (A2HDM) pQ, where each 
pair of Yukawa coupling matrices (describing the couplings to right-handed up-type quarks, 
down-type quarks and leptons, respectively) is assumed to be aligned in flavour space. The 
resulting Yukawa Lagrangian [T], 

C Y = H + (x) [ft VM d V R - MlVV L ] d(x) + q v{x)M{P R l{x)} 



~ E vfviW f(x)M f V R f(x) + h.c. , 



v 
<p,f 

contains three complex couplings $f (/ = u, d, I) which parametrize all possible freedom 
allowed by the alignment conditions. Vr,l — 1 ^ 75 are the right-handed and left-handed 
projectors and the couplings of the three neutral scalar fields ${x) = {h(x), H(x), A(x)} are 
given in terms of and the orthogonal transformation 1Z relating the scalar weak and mass 

eigenstates: = Tin + (H i2 + iTZ i3 ) ft,z, Uu = Kil + \R>i2 ~ if^m) <£■ 

While the CKM matrix V remains the only source of flavour-changing phenomena, the 
universal (flavour-blind) couplings introduce three new complex phases and, therefore, 
a new source of CP violation. For particular (real) values of these parameters the usual 
CP-conserving models based on discrete Z-i symmetries are recovered. Quantum corrections 
induce a misalignment of the Yukawa matrices, generating small FCNC effects pT|2|20|21 j . The 
flavour symmetries of the A2HDM strongly constrain the allowed FCNC structures, providing 
at the quantum level [T][2] an explicit implementation of the popular Minimal Flavour Violation 
(MFV) scenarios [22-429] . but allowing at the same time for new CP-violating phases. 

A MFV structure within the context of the type II two-Higgs doublet model (2HDM) was 
considered in |25] and flavour-blind phases have been recently incorporated in [28jl]. These 
references perform a perturbative expansion around the usual U(1)pq symmetry limit (type 
II) and look for tan /3-enhanced effects. Since the A2HDM does not assume any starting 
ad-hoc symmetry, it leads to a more general MFV framework with tan /3 substituted by the 
six-dimensional parameter space spanned by the couplings $f. While giving rise to a much 
richer phenomenology, the A2HDM implies an interesting hierarchy of FCNC effects, avoiding 
the stringent experimental constraints for light-quark systems and allowing at the same time 
for potentially relevant signals in heavy-quark transitions [2J. Notice that in the general case, 
without U(1)pq or Z 2 symmetries, tan/3 does not have any physical meaning because it can 



Regarding MFV and CP violating phases, see also 27,30 -33 . 
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be changed at will through SU(2) field redefinitions in the scalar space; the physics needs to 
be described through the (scalar-basis independent) parameters 

The presence of new weak phases immediately poses the question of compatibility with 
existing measurements of CP violating observables. One decay where a high sensitivity is 
expected is b — >■ s~f, for three reasons: the SM asymmetry is known to be tiny [34H36] . 
its measurement is rather precise and compatible with zero, and the potential influence of 
new physics in this decay is large. In addition, being an inclusive process, non-perturbative 
uncertainties are expected to be under better control than for exclusive modes. The branching 
ratio and the CP rate asymmetry have both been analyzed within the context of a general 
2HDM [37], and it has been pointed out that in large regions of the parameter space the 
NLO predictions are not reliable, because they suffer from a very large renormalization- 
scale dependence. For the part of the parameter space not affected by these problems the 
asymmetry was found to be small. 

In this work, we update the calculation for the asymmetry within the framework of 
the A2HDM, extending the phenomenological analysis and discussing its implications for 
the model parameters. We will show that the combined phenomenological constraints from 
Z — > bb, S-mixing, -mixing and Br(6 — >■ 57) [2j exclude the problematic regions of the pa- 
rameter space. The predicted range for the CP asymmetry turns out to be within the present 
experimental limits for all allowed ranges of the A2HDM couplings, although for some par- 
ticular values of the parameters it is not far from the achieved sensitivity. This makes future 
high-precision measurements of this quantity very desirable and calls for a better control of 
theoretical uncertainties. 

To that end, we first analyze the branching ratio in some more detail. We have started 
to discuss this observable in [2j together with others, but deem a more detailed analysis in 
this context worthwhile, due to its relevance for the rate asymmetry, but also due to its own 
high sensitivity and its intimate relation to the like-sign dimuon charge asymmetry (LDCA) 
in Z?°-decays. 

We proceed as follows: In section [2] we discuss the B —> X s ^f branching-ratio prediction 
within the A2HDM, extending our results derived in [2]. In section [3] we analyze the CP rate 
asymmetry, with emphasis on its scale dependence, and discuss its potential impact on the 
parameter space of the A2HDM. Finally, we conclude in section [4] with a brief summary of our 
findings, and analyze their impact regarding the possible influence of charged-scalar effects 
on the LDCA in our framework. 

2 B — > X s 7 branching ratio 

Because of its importance for the following discussion, we first perform an analysis of the 
B —7- X s "f branching ratio. In the SM, this decay is known to be subject to large radiative 
corrections. While the problem of a sizable scale dependence has been basically resolved 
with the calculation of the NLO corrections, the issue of charm-mass scheme dependence can 
only be addressed at NNLO [38]. To achieve this, a huge effort is being made, and by now 



2 



the branching ratio is essentialljQ calculated up to NNLO in the SM. We follow here the 
calculation by [42], giving Br(B — > X s j) E ^ g° GeV = (3.15 ± 0.23) x 10~ 4 , in agreement with 
the present world average [43] 

Br (5 -> X s7 )f 6GeV = (3.55 ± 0.26) x 10~ 4 . (2) 

Note however, that different treatments of photon-energy-cut related effects lead to slightly 
different results [44H46] (see also [47]). The related shifts are of the order of the uncertainty 
assigned in [42]. Regarding the non-perturbative part of this calculation, contributions with 
the photon coupling to light partons ( "resolved" photon contributions) lead to the appearance 
of non-local matrix elements, implying an unreducible error of ~ 5% [48] . 

In the 2HDM, the calculation has been performed up to NLO [37U49ff51~] . For m s = 0, 
the effective low-energy operator basis remains the same as in the SM and the modifications 
induced by new-physics contributions appear only in the Wilson coefficients: 

Cf (nw) = Ci t sM + k«| 2 Ci )UU — Ci tU d , (3) 

where <f*^ = \s u \ \$d\ ^ defines the relative phase ip. Values around ip = (tt) correspond 
to destructive (constructive) interference with the SM amplitude, as occurs in the 2HDM of 
type I (II). 

Since the 2HDM contribution to the decay amplitude has only been calculated up to 
NLO, the corresponding terms of 0(a^) are neglected consistently in the rate. An important 
observation in [37] is that the stabilization against scale variations does not work as well for 
the general 2HDM as for the SM, because in some regions of the parameter space cancellations 
between SM and new-physics contributions occur, enhancing the sensitivity to higher-order 
QCD corrections. In [2] we have used the NNLO SM calculation, retaining terms up to order 
a 1 sm an d a s,NP which stabilize the result, however with some problematic regions, leading 
to unphysical results, still remaining. But these regions (corresponding to large values of |? u |) 
are excluded by constraints coming from Z — > bb, ex and Am B o, allowing us to perform a 
consistent theoretical study of Br(.B — > X s ^f) in the whole remaining parameter space. 

In Fig. [I] we show the resulting constraint in the plane — |q|, varying the remaining 
parameters in the ranges M H ± 6 [80,500] GeV and if G [0, 2ir], and treating the errors as 
described in [2j. Due to the additional degrees of freedom, the resulting constraint is relatively 
weak; especially, no direct bound can be found with respect to the charged-scalar mass, in 
striking contrast to the 2HDM type II. As discussed in [2], this changes as soon as the phase 
is kept fixed to a certain value. 

Numerically, the decay amplitude has roughly the following structure for large scalar 
masses (using numerical estimates for the different contributions from [37]): 

a , (\ n , * /500 GeV\ 2 nM . l2 /500 GeV\ 2 l 
A^M|l-0.1^^--j +0.01 ) • (4) 

2 For the present status and recent developments, see e.g. pMUT] . 
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Figure 1: The constraint from Br(6 — > s'j) plotted in the plane \<^ u \ — |^|, see text. 

From this it is obvious that we can expect in that case constraints on the parameter combi- 
nations |<f u | 2 /Mo- ± and ?*<^/Mjl ± , the latter being complex. For sizable |<^|, the last term 
is negligible, as is constrained to be 0(1) at most, leaving only the dependence on the 
combination <r*c^/M^ ± . The resulting limits on the single parameters in this ratio are rel- 
atively weak, as seen above. The strength of the constraint lies however in creating strong 
correlations between different parameter combinations. This is illustrated in Fig. [21 where 
the product \^d\ is plotted against the charged-scalar mass M H ± and the relative phase (p, 
respectively. The hole in the left plot can be understood as separating two regions where 




1 2 3 4 5 6 100 200 300 400 500 
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Figure 2: The constraint from Br(6 —¥ sj) plotted in the planes V 3 O e ft) an< ^ \'*u$d\~Mff± 

(right). 

the NP influence is relatively small (lower part) and where it is approximately twice the size 
of the SM contribution (upper part), the latter corresponding to a fine-tuned solution. For 
larger phases it is not possible anymore to cancel the SM amplitude completely, so only one 
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constraint is remaining, making the gap disappear. From the right plot it becomes obvious 
that high values for the product are only allowed for large values of the charged-scalar 

mass, due to the restriction on the ratio. The effect of these correlations can obviously be 
large; it is taken into account in the following analysis. 

3 CP rate asymmetry 

The CP rate asymmetry for the process b — > S7 is defined as 

_ Br(£ -» X sl ) - Br(B -» Xfj) 
aCP Br(B -> X sl ) + Br(B -> X- sl ) ' [ > 

Being doubly Cabibbo suppressed, it is known to be tiny in the SM, below 1% [34ti36] . making 
it a sensitive probe of new-physics effects. The relative scale uncertainty is large, as it is to 
be expected from a NLO calculation for a CP asymmetry, because it is the first contributing 
order. However, it is small in absolute terms. The same is true for the charm-mass dependence. 
Note that the enhanced power-corrections at order 0(a s ^-) [38], mentioned above in the 
discussion of the rate, are not yet included in the calculation for the asymmetry. Being of 
relative order A/m&, they are likely to increase the uncertainty significantly. However, the 
experimental uncertainty of the present world average |43j . 

Oc|(B^l s 7) = -0.012 ±0.028, (6) 

is much larger than the expected value, so the calculation at NNLO has been considered less 
interesting up to now. 

Within the context of the general 2HDM the analysis has been performed at NLO [37J, 
working in the limit V ub V* s = 0, which implies by unitarity that V cb V* s = — V tb Vt s . It was 
pointed out in [37] that, given the cancellation problems in the branching ratio described 
before, the prediction was not reliable in part of the parameter space where it exhibited a 
large scale dependence, but that for parameter choices where the branching ratio was well 
behaved the predicted acp was very small, 0(1%). 

We now turn to reanalyze the asymmetry in the context of the A2HDM. In our previous 
study we have used the SM NNLO calculation of the branching ratio; however, we will not 
include NNLO information in the predicted asymmetry, for two reasons: First, the necessary 
calculation does not exist completely, and the interpolation for the charm-mass dependence 
has been done specifically for the branching ratio, ignoring some imaginary parts which can 
be relevant for the asymmetry. Second, and more importantly, in our approximation the 
asymmetry is a pure new-physics effect. Therefore, QCD corrections should enter at the same 
level for the SM and new-physics amplitudes. Moreover, as the scale dependence at NLO is in 
general more severe for the A2HDM than for the SM, inclusion of the NNLO SM corrections 
only would not stabilize the full result. 

We confirmed these expectations explicitly by calculating the asymmetry with the NNLO 
SM contributions to the branching ratio included. The result showed a relatively large shift 



5 



in the central value, which we attribute to the different charm-mass dependence, and no sta- 
bilization of the scale dependence at all. This behaviour is uniform for different values of q U) d 
and Mjj± . We conclude that for the inclusion of NNLO corrections important parts are miss- 
ing, the calculation of which is beyond the scope of this paper. This lack of NNLO corrections 
will make the problems described in the branching-ratio analysis reappear. Nevertheless, we 
start by analyzing the asymmetry for the same range of parameters considered in our previous 
study of the branching ratio, |<j u | G [0,2], \q d \ £ [0,50], <p G [0,2n], and M H ± e [80,500] GeV, 
to examine the strength of this observable by itself. 

Figure [3] shows the [i^ renormalization-scale dependence of the predicted CP asymmetry 
at NLO, as a function of the phase <p, taking M H ± = 200 GeV, \s u \ = 0.1 and \$d\ = 5 
(left) and 50 (right). The black central curve represents acp for m, = 2.5 GeV, while the 
outer (larger absolute values) and inner (smaller absolute values) lines correspond to fit, = 2 
and 5 GeV, respectively, using the same range of variation considered in the branching ratio 
analysis. We observe that the //& dependence is proportional to |acp|; therefore in the regions 
where the asymmetry is relevant the theoretical error from the scale dependence is as well. 
As noted before, the overall scale dependence is strong, approximately 25% in most of the 
parameter space. If not noted otherwise, in the following the scale is fixed to = 2.5 GeV. 

4=0.1, Q=5, M H =200 GeV £,=0.1, 4=50, M w =200 GeV 




(p/md (p/md 



Figure 3: CP asymmetry at NLO, as function of the relative phase tp, for q u = 0.1, q = 5 
(left) and 50 (right) and M H ± = 200 GeV. The band shows the variation with the scale, taking 
fib = 2.5 GeV (central black line), 2 GeV (outer line) and 5 GeV (inner line). 

The asymmetry is small for most values of the phase ip, apart from values around zero, 
where for some parameter choices one observes an unphysical enhancement of |acp|- This 
is apparent in Fig. [3] (right panel), corresponding to q = 50, which shows a large predicted 
asymmetry around ip = ±0.2 rad (modulo 2n), with a very large scale dependence. This 
enhancement does not correspond to any large CP contribution, but rather to a destructive 
interference in the decay amplitude leading to a nearly vanishing rate, which is only possible 
for small imaginary parts. This resonant behaviour can again be easily understood from Eq.[H 
a zero of the amplitude occurs around ?*?d/Mj^ ± ~ 4 x 10 -5 GeV -2 with e lip ~ 1. For lower 
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scalar masses the \s u \ 2 term can give a more sizable contribution, modifying this simplified 
behaviour, but again the amplitude can become small through the cancellation of the SM 
and new-physics contributions. We illustrate this fact in Fig. HI which shows the parameter 
region in the |<j u ||^|/Mjl ± — M H ± plane leading to values of acp outside the (95% CL) allowed 
experimental range, taking ip = 6.15 rad and fixing the renormalization scale at \x\, = 2.5 GeV. 

tp = 6. 15 rad 

0.0006 
0.0004 
0.0002 

100 200 300 400 500°' 

M H I GeV 

Figure 4: Region in the |? U ||^|/M^ ± — M H ± plane leading to values of acp outside the (95% 
CL) allowed experimental range, for fixed phase ip = 6.15 rad and scale fib = 2.5 GeV (M H ± 
in GeV units). 

The regions of small amplitudes around values of p ~ 0, 2tt are of course very sensitive to 
the adopted truncation of the perturbative expansion in powers of the strong coupling and 
introduce a correspondingly large theoretical uncertainty. Their appearance signals the need 
to incorporate higher-order corrections into the calculation. Fortunately, since the measured 
branching ratio is definitely non zero (by many standard deviations) and compatible with the 
SM prediction, these problematic regions of the parameter space, where cancellations occur, 
are already excluded. Imposing the constraint that the experimental branching ratio should 
be reproduced, the surviving allowed parameter ranges lead to well-behaved amplitudes and 
small asymmetries below the present limits. 

In Fig. [5]we plot the predicted maximal value for the CP rate asymmetry, as a function 
of the Yukawa phase p = arg(?*? rf ), including the Br (B — > X s ^f) constraint by showing only 
points that lie within its 95% CL experimental range. We have scanned the predictions over 
the parameter ranges M H ± G [80,500] GeV, |? M | G [0,2] and |<^| £ [0,50], and have treated 
the error for the branching ratio as explained in [2] . We find as expected that the problematic 
regions have been excluded, and that the maximal achievable asymmetry is compatible with 
the present experimental limits at 95% CL, given the scale dependence of the prediction. 
Sizable asymmetries at the 1-5% level seem possible for p ~ ±0.7 rad, but the corresponding 
theoretical uncertainty is unfortunately quite large, as shown by the large scale-dependence 
of the theoretical results. 
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Figure 5: Maximal value of acp versus the phase ip at NLO, for Mjj± E [80,500] GeV, 
E [0,2] and |q| E [0,50], taking into account the experimental Br(B — > X s j) constraint at 
95% CL. The three curves correspond to fib = 2 (outer), 2.5 (center) and 5 GeV (inner). The 
dotted (continuous) horizontal lines denote the allowed experimental range at 95% CL (68% 
CL). The vertical shadowed band indicates the region in which an enhancement of the LDCA 
induced by charged-scalar effects in the direction of the measured value may occur. 

4 Discussion 

The A2HDM introduces new sources of CP violation in the flavour sector while avoiding 
FCNCs at tree level, and provides an explicit counter-example to the widespread assumption 
that in 2HDMs without tree-level-FCNCs all CP- violating phenomena should originate in the 
CKM matrix. Since all Yukawa couplings are proportional to fermion masses, the A2HDM 
gives rise to an interesting hierarchy of FCNC effects, avoiding the stringent experimental 
constraints for light-quark systems and allowing at the same time for interesting signals in 
heavy-quark transitions. The flavour-blind phases present in the model imply a broad range 
of very interesting phenomenological implications, opening new possibilities which are worth 
to be investigated. 

An obvious question to address is the phenomenologically allowed size for the new phases, 
which should be constrained by the measured CP-violating observables in flavour physics. 
Of these, the b — > sj asymmetry is known to have both a very small SM contribution and 
a high sensitivity to new physics. Unfortunately, although the b — > sj rate is known at 
NNLO in the SM, the 2HDM contributions have been only computed at NLO. Since a non- 
zero rate asymmetry requires absorptive contributions leading to a strong-phase difference, 
it first appears at NLO. Therefore, the theoretical uncertainties are much larger than for 
the branching ratio. On the other hand, some of the parametric uncertainties cancel in the 
asymmetry. 

We have analyzed the predicted b — > sj rate asymmetry within the A2HDM, scanning the 
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parameter space of the model over the full domain allowed by known constraints from other 
processes, which we have derived in a previous work [2]. We have shown that, while acp 
could be enhanced for some particular values of the parameters, it remains compatible with 
its present measurement when considering parameters leading to acceptable values for the 
branching ratio. Sizable CP asymmetries at the 1-5% level, close to the present experimen- 
tal bound, seem possible for Yukawa phases around arg(?*<^) ~ ±0.7 rad, making a future 
measurement of this quantity very interesting. Such a measurement would be possible at a 
Super- B factory, where precisions better than 1% could be achieved [52j[53] . However, the 
presently large theoretical uncertainty from scale and mass dependence makes necessary a 
complete NNLO calculation for the SM and the A2HDM contributions to fully exploit such 
a measurement. At the 1% level, one also needs to analyze the small corrections induced by 
|T4feKTsl 7^ 0) which have been neglected up to now. An exception is given by potential exper- 
imental information on the quadrant of the new-physics weak phase (p: once a measurement 
is available which determines the sign of the asymmetry unambiguously, two quadrants are 
excluded. 

We also extended the phenomenological analysis of the branching ratio over the first as- 
sessment done in our previous publication [2] . The strong constraints on the model parameters 
e.g. in type II models are not visible in the A2HDM, due to the presence of additional param- 
eters. They are replaced by correspondingly strong correlations, which we analyzed in some 
detail above, and which play an important role when combining other constraints with this 
observable. 

An example for their influence is given by the LDCA in £>°-decays, which with the aid 
of the corresponding measurement at the i?-factories for B® only can be interpreted as an 
enhanced asymmetry in the -B^-system, 

a s sl = lm[ T ^-)= sin</> s = tan0 s , (7) 



Mf 2 J \Mf 2 \ AM B o 

where M? 2 — | Tf 2 = (B®\H^ S B = 2 \B°) and (j) s = arg(-Mf 2 /rf 2 ). Noting that in the A2HDM 
the change in the rate difference is negligible, we can parametrize the relative influence of the 
charged scalar on the asymmetry as 

where Mf 2 |f u n = Mf 2 |sMA s . The possible range for this ratio is shown in Fig. [6] (details on 
the calculation can be found in [2]): the plot on the left hand side takes into account the 
constraint shown in Fig. [TJ varying all NP parameters in the ranges given before. However, 
taking the correlations into account by using a parametrization |<j*Q| max (y?, M^±) for the 
constraints shown in Fig. [21 the potential influence of charged-scalar effects on this observable 
is reduced by roughly a factor of 10 compared to the result without correlations. While an 
enhancement by a factor of ~ 5 is still possible, larger effects from charged-scalar contributions 
are excluded. This large impact can be understood as due to a combination of the two 
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Figure 6: Possible influence of charged- scalar- effects on the LDCA, taking the constraint from 
Br(6 —7- 57) shown in Fig. [7] (left) and Figs\^ (right) into account. Shown is the allowed range 
for the relative factor with respect to the SM result as a function of the phase ip; see for 
details of the calculation. 

correlations shown above: as discussed in [2 J , sizable values of the LDCA require relatively low 
charged-scalar masses, a sizable phase and large values for the product |<^?d|- The correlations 
in the constraint from 6—7-57 imply that the latter can only be achieved for relatively 
high masses and small values of the phase. Therefore the full suppression shows up only 
in the correlated analysis performed here. A possible enhancement from additional neutral 
scalars remains untouched by this. The present data suggest a preferred negative sign for the 
B® semileptonic asymmetry a s sl , which for charged-scalar effects to be in the right direction 
requires arg(?*<^) E [37r/2, 27r]. Due to the different structure of its amplitude, the b — > S7 
rate asymmetry could provide complementary information on the Yukawa phase. 

Note added 

After this paper was posted in the arXiv, a very relevant work on the rate asymmetry has 
appeared [S3], showing that it receives a (parametrically leading in the SM) long-distance 
contribution arising from the interference of the electromagnetic dipole amplitude with an 
up-quark penguin transition accompanied by soft gluon emission. 
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